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Abstract 

The motion of ions in the fields B = Bqz and E = E^y cos(k±y — u>t) is con- 
sidered. When u ^> fli and v± > uj/kj_, the equations of motion may be reduced 
to a set of difference equations. These equations exhibit stochastic behavior when 
Eq exceeds a threshold. The diffusion coefficient above the threshold is deter- 
mined. Far above the threshold, ion Landau damping is recovered. Extension of 
the method to include parallel propagation is outlined. 
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Equations of Motion 

Consider an ion in a uniform magnetic field and a perpendicularly propagat- 
ing electrostatic wave, 

B = B a z, E = E y cos(k±y - uit). (f) 

Normalizing lengths to fc^ 1 and times to fi" 1 (f2, = qiBo/rrii), the Lorentz force 
law for the ion becomes 

y + y = acos(y - vt), x = y, (2) 

where 

v = uj/fli, (3a) 

a= w (3b) 

We solve (2) by approximating the force due to the wave by impulses at 
those points where the phase is slowly varying, i.e., at y = v. The trajectory of 
the ion is given in Fig. 1. Expanding the trajectory about the resonance point, 
we find that the magnitude of the impulses is given by 



/oo 
acos(0 c - \y c t' 2 )dt' 
-CO 



= a{2ir/ \y c \) 1/2 cos[<?!> c - sign(y c )7r/4], (4) 

where 4> c = y c — vt c and t c and y c are the time and position of the wave-particle 
"collision." We may determine the Larmor radius and phase of the ion at the end 
of the jth orbit [the beginning of the (j + l)th orbit] in terms of these quantities 
at the beginning of the jth orbit. (Details are given in Ref. 1.) The resulting 
difference equations are 

u = 9 — p, v = 8 + p, (5a) 

6 = \{v + u), p = \{v-u), (5b) 

Uj + i — uj = 2nd — 2irAcosvj, (5c) 

Vj+i — Vj = 2nS + 2ttA cos u^+i , (5d) 



where 



9 = ut (mod27r), (6a) 

p = (r 2 - v 2 ) 1 ' 2 - !/cos _1 (z//r) + vtt - tt/4, (6b) 

A- $"2*^. (7a, 

8 = v-n. (7b) 
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Here r is the normalized Larmor radius, k±v±/£li, and n is an integer. The limits 
of validity of (5) are 

r-^>(i^) 1 / 3 , A« (r 2 -v 2 f' 2 /r 2 . (8) 

In Fig. 2 we compare the trajectories obtained using the exact equations of 
motion, (2), with those obtained from the difference equations, (5). We see that 
the agreement is very good indicating that (5) is an excellent approximation of 
(2). 

There are three advantages to using the difference equations in preference 
to the Lorcntz force law. Firstly, they are much quicker to solve numerically. 
Secondly, because of the way the equations were derived, the results are easier 
to interpret. Lastly, the equations have separated out two velocity-space scales, 
the p scale (~ Qi/k±) and the r scale (~ w/k±). We therefore treat A which 
is a function of r as a constant when iterating the equations. This means that 
the diffusion coefficient is independent of p and so is much easier to determine 
numerically 



Examples of Trajectories 

When A is infinitesimal, (5) may be solved by integrating (summing?) over 
unperturbed orbits. Substituting Uj = uq + 2it8j and vj = vo + 2ir8j into the 
right hand sides of (5c) and (5d) gives 

Pn — po = 2irAcos(po — n5) 



x < 



(—1) cos9 N, for S = integer, 

sin(27r<W + O ) - sin O , r , . ( 9 ) 

, for 5 ^ integer. 

2 sm(7rd) 

Note that the trajectory is secular or not depending on whether or not S is 
an integer (uj is a cyclotron harmonic). Formally, we may compute a diffusion 
coefficient using 

V= Lim «^-^) 2 >. (10) 
Substituting (9) into (10) gives 

oo 

T> = ir 2 A 2 cos 2 po S(S-m), (11) 

m— — oo 

where 5 is the Dirac delta function. Converting back to r and t and undoing the 
normalizations gives 

/ m = — oo x ' 

which is the usual quasi-linear diffusion coefficient. 

If we consider finite but small A, then all the trajectories are bounded. There 
are three distinct cases, S — (which is the case considered by Fukuyama et al. 2 ), 
5 = |, and 5 ^ or |. The trajectories for these cases are shown in Fig. 3. 
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When A is increased, the system undergoes a stochastic transition, an ex- 
ample of which is shown in Fig. 4 for S = 0.23. Below the stochasticity threshold, 
nearly all the trajectories are integrable [Fig. 4(a)] or, if there are stochastic tra- 
jectories, they are bounded in p [Fig. 4(b)]. Above the threshold, nearly all the 
trajectories are stochastic and unbounded. The value of the threshold may be 
numerically determined and is found to be A = A s = j . Above this value of A, 
the kick received by the ion during one transit through resonance is sufficient to 
change the phase of the kick received when next in resonance by ir/2. 



Diffusion Coefficient 

When computing the diffusion coefficient numerically, it is convenient to 
work with the correlation function, Ck, where 

C k = (cija-j+k), (13) 

cij is the particle acceleration, aj — Pj+i — Pj, and the average is over an ensemble 
of particles and over the length of a given trajectory (i.e., over j). Then the 
diffusion coefficient, V, is given by 

oo 

V=-C + J2c k . (14) 
1 k=i 

[This definition is equivalent to (10).] The advantages of defining T> in this way 
are twofold. Firstly, the statistical fluctuations in the computation are minimized. 
Secondly, it is easy to introduce the effects of collisions on the diffusion coefficient. 
This is accomplished as follows: If fco is the mean number of cyclotron periods 
between decorrelating collisions (such collisions need only result in deflection by 
l/v, a small angle), then the probability of such a collision taking place in k 
periods is 1 — cxp(— k/k ) since collisions are independent events. Collisions 
may then be included in the computations of V by replacing Ck in (14) by 
C k exp(—k/ko). 

In the limit A ^> A s , the kicks the ion receives are uncorrelated so that only 
Co is nonzero. Assuming that the trajectory is ergodic, we obtain T> = ^ir 2 A 2 . 
When A is not large, we account for the correlations between the kicks received 
by the ion by writing 

V=\n 2 A 2 g 2 {A), (15) 
Numerically determining g(A) we find that approximately 

g(A) =m^(l-A 2 JA 2 ,0), (16) 

with A s = \. 

Converting (15) back into usual variables we obtain 

-_£]%( u V g\- 

2 ml \k±v±) (k\vl -w 2 )V2- 

In the limit A ^> A s , when g(A) w 1, this is just the zero-magnetic-field result, 
{-K/2){qi/mi) 2 E 2 5{lo — k±v y ), averaged over Larmor phase. Thus, in this limit, 
we recover ion Landau damping. 



D= im(jjL.y 9{A) d7) 

u o ™2 (Li „.2 ..2M/2- y 1 ') 
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Extension to Parallel Propagation 

The diffusion coefficient was so easily calculated above because of the simp- 
lification obtained by reducing the problem to difference equations, (5). This 
reduction may be achieved in similar problems. We consider here the case where 
the wave has some component of parallel propagation so that (1) becomes 

B = B z, E = Eo(y + k\\z/k±) cos(k±y + k\\z - ut). (18) 

Adopting the same normalization as before, we obtain 

y + y = acos(y + (z — ut), (19a) 

z = a£cos(y + C,z — vt), (19b) 

where ( — k\\/k±. The difference equations for a particle with normalized Larmor 
radius, r, and normalized parallel velocity, w — are 

u = 6-p, v =9 + p, (20a) 

Uj + i — Uj = 2n-fj +1 / 2 — 2-irAcosvj, (20b) 

Vj+i — Vj = 2n-fj + i/ 2 + 2wAcosUj + i, (20c) 

Tj+i/2 = S- 0{pj + irAcosvj), (20d) 
where the variables and p are given by 

9 = vt-(z (mod 2?r), (21a) 

p=(r 2 - /i 2 ) 1/2 - Atcos _1 (/i/r) + pir - tt/4 - 2-Km. (21b) 
The definitions of the parameters A, j3, and p are 



2\ 1/2 apQ 



(22a) 



nj r(r 2 -^ 2 )V4' 

/3 - C 2 r/(pQ), (22b) 
S = p + fip-n. (22c) 
Here, m and n are integers and 

p = v — (w, (23a) 



Q = 



(r 2 -p 2 ) 1/2 f _!/Ml( 2 r 



7T — COS 



- (?)] 't- (m) 

Despite appearances 5 is a parameter independent of p since the quantity p + /3p 
is a constant. (This follows from energy conservation in the wave frame.) The 
restrictions on the validity of (20) are 

r- p->{\p) 1/z . (24) 
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A comparison between the exact equations of motion, (19), and the difference 
equations, (20), is shown in Fig. 5 for ( = 1 (propagation at 45°). Again, there 
is excellent agreement. 

The results of Smith and Kaufman 3 may be obtained in the limit (3 —* oo 
and A — > (Q — > 0). In that case, the change in p is negligible so that it is 
necessary to rescale the velocity variable by defining 

o j = Airfipj - 2tt<5. (25) 

Equation (20) then becomes 

6 j+1 - 6j = -<7j - \Kcos{6j + p), (26a) 

CTj+i - (Tj = \K[cos(9 j+ i - p) + cos(9j + p)], (26b) 
where K = 4:ir 2 Af3. In (26) p is a constant. Setting 

i> = a + \Kcos(d + p) (27) 

gives 

e j+1 -e j = -^, (28a) 

tpj+i — tpj = K cos p cos Oj+i . (28b) 

This is the "standard mapping" studied by Chirikov. 4 The island overlap condi- 
tion for this mapping is |_ftTcosp| > 7r 2 /4 or 

\a\ > |l6C 2 J M (r)| _1 , (29) 

which is the stochasticity threshold obtained by Smith and Kaufman. The island 
overlap criterion is a significant overestimate of the stochasticity threshold for the 
standard mapping. 4 Greene 5 has calculated that the true threshold is a factor 
(7r 2 /4)/0. 971635 w l\ smaller than the result given above. 
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Fig. 1. Motion of an ion in velocity space, showing the kicks it receives when 
passing through wave-particle resonance. 
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Fig. 2. Comparison of the difference equations with the Lorentz force law. (a) 
Trajectories computed using (2) with v = 30.23 and a = 2.2. (b) Trajectories 
computed using (5) with 5 = 0.23 and A — 0.1424, which are given by (7) with 
v = 30.23, a = 2.2, and r = 47.5. In each case the trajectories of 24 particles are 
followed for 300 orbits. 




Fig. 4. The trajectories of ions with 5 = 0.23 and (a) A = 0.05, (b) A = 0.2, (c) 
A = 0.35. The initial positions are shown by crosses. 
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Fig. 5. Comparison of the difference equations with the Lorentz force law for 
finite fen. (a) Trajectories computed using (19) with v = 20.23, a = 0.23, and 
( = 1. The total energy (kinetic + electrostatic) of each particle is the same 
and is chosen so that w ~ when r = 31. (b) Trajectories computed using 
(20) with 5 = 0.8548, A = -0.06768, and (3 = -0.5595 which are given by (22) 
with v = 20.23, a — 0.23, £ = 1, r = 31, w = 0, and m = 10. In each case 
the trajectories of 24 particles are followed for 300 orbits. The vertical scale is 
inverted in (b) since Q < 0. 



